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Overview 

For a few fleeting moments after the Big Bang the universe was filled with an astonishingly hot and dense 
soup known as the Quark-Gluon Plasma (QGP), a precursor to the matter we observe today that consisted of 
elementary particles. Although this Quark-Gluon Plasma also contained leptons and weak gauge bosons, its transport 
properties were dominated by the strong interaction between quarks and gluons. Utilizing the most powerful particle 
accelerators, physicists now conduct head-on collisions between heavy ions, such as gold or lead nuclei, to recreate 
conditions that existed at the birth of the universe. The most striking discovery in relativistic heavy-ion collisions is 
that the hot and dense matter created during the collisions behaves like an almost perfect (inviscid) liquid, meaning 
that it can be characterized by a very small shear viscosity (77) to entropy density (s) ratio (or, “specific shear 
viscosity” )m- It turns out that p/s for the QGP is smaller than that of any known substance, including that of 
superfiuid liquid helium. In Fig. we illustrate schematically the specific shear viscosity p/s normalized by 
the minimum bound in a large class of theories with infinitely strong coupling [ 5 ], for four different types of fluids. 
The QGP at high temperature exhibits the smallest value of 77/5 of any fluids occurring in nature. 

Since the discovery around the turn of the millenium of the QGP and its surprisingly perfect fluidity, strong 
interest has emerged in both theoretical and experimental work to constrain the transport properties of the QGP 
in relativistic heavy-ion collisions. The specific shear viscosity characterizes one of the most important transport 
properties of the QGP. It is at present impossible to compute this transport coefficient with any sort of precision 
from first principles. Alternatively, theoretical analysis of the experiments conducted at the Relativistic Heavy- 
Ion Gollider (RHIG) and the Large Hadron Gollider (LHG) offers opportunities to unravel phenomenologically the 
collective behavior of the QGP created at high energies. Importantly, since the fireballs created in these “Little 
Bangs” are minuscule {V ^ and cool almost instantly 5 x physicists can only study the QGP 

through its remnants. This raises unique experimental and theoretical challenges. To trace the information of the 
detected stable particles back to the QGP phase at the early stages of the collisions, in order to extract the QGP’s 
properties, requires a quantitative understanding of the entire evolution of a heavy ion collision, from the formation 
and thermalization of the QGP to the dynamics of the hadron resonance gas into which it eventually decays. [A 
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FIG. 1. The fluid imperfection index of various fluids as a function of temperature. This picture is taken from Tribble 

R (chair), Burrows A et al. 2013 Implementing the 2007Long Range Plan, Report to the Nuclear Science Advisory Committee, 
January 31, 2013. Available at http://scieiice.energy.gov/np/nsac/reports/. 
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FIG. 2. Left Panel: Illustration of the dynamical evolution of relativistic heavy-ion collisions. Right Panel: The collective 
space time evolution of the collision systems created at the RHIC and the LHC compared with ultra-cold fermi gas [9]- The 
temperature is color coded. 


schematic view of the evolution as we currently understand it is illustrated in the left panel of Fig.|^] Such a reverse 
engineering process requires good control of all the model uncertainties, in addition to a consistent description of a 
large variety of measured observables. 

In this featured article we highlight the major recent advances in theory and experiment in constraining the 
QGP shear viscosity. We begin by explaining the main idea how to measure rj/s in heavy-ion collisions. 

The QGP viscometer 

• Bulk dynamics and charged hadron elliptic flow: the state of the art 

The success of viscous hydrodynamics in describing the bulk evolution of relativistic heavy-ion collisions at 
RHIC [6] and predicting and describing flow data later collected at the LHC miH] leads to an important conclusion: 
At the macroscopic level, initial anisotropic pressure gradients in the fireball drive the system to collectively develop 
a momentum anisotropy, as illustrated in the right panel of Fig. This dynamic evolution shares a great deal of 
similarity with other strongly coupled many-body systems, such as ultra-cold fermi gases [ 9 ]. Inhomogeneities and 
anisotropies in the initial density distribution can be characterized by the spatial eccentricities ^n} defined by 

. Jrdrd(l>r^e{r,(t>)e^<l’ . Jrdrd<Pr^e{r, ^ 

^ J rdrd(pr^e{r,(f>) " ^ rdrd(j)r'^e{r,(j)) ~ ’ 


where r and (j) are polar coordinates in the transverse plane perpendicular to the beam direction. (At top RHIC and 
LHC energies, the dynamics along the beam direction is, to good approximation, boost-invariant, i.e. independent 
of the longitudinal motion of the reference frame, and therefore less informative for the present discussion than the 
transverse dynamics.) Similarly, the azimuthal dependence of the emitted particles’ momentum distribution can be 
characterized by the anisotropic flow coefficients Vn and their associated event plane angles defined by 


Vn 




U ' 


J d(pp 


and 


f d(j)p 


dN 

dypTdpTdcf)p 


^incpp 




dN 


P dypTdpTd(j)p 


( 2 ) 


where pr and <f>p are polar coordinates in transverse momentum space. The conversion efficiency from the initial 
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FIG. 3. Eccentricity-scaled charged hadron elliptic flow as a function of final charged hadron multiplicity per overlapping area 
compared to a viscous hydrodynamics + hadron cascade hybrid simulation in Au+Au collisions at top RHIC energy. The plot 
is taken from m- 


{sn} to the final {t>„} is degraded by viscosity in the medium m-, this is the main effect that opens the door to 
measuring r]/s with anisotropic flow observables. 

In particular, the second order anisotropic flow coefficient, dubbed the elliptic flow V2^ of charged hadrons was 
first measured in heavy-ion experiments at the CERN SPS and RHIC. In Fig.[^ we compare the RHIC measurements 
to results from viscous hydrodynamic simulations, in order to extract 77/s of the QGP. By confronting state-of-the-art 
phenomenological modeling with precise experimental measurements, the specific shear viscosity of the QGP can be 
constrained as m 
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A glance at Fig.[ 2 shows that even the upper limit lies much below the minimal specific shear viscosities seen observed 
in any other previously known fluid. However, since the central value is so small, the remaining relative uncertainty 
is still large. It is dominated by model uncertainties in the initial conditions (the two panels in Fig. show the same 
data compared with two typical initial-state models); it is irreducible unless further anisotropic flow observables, in 
addition to the elliptic flow, are taken into account. 

• Event-by-event anisotropic flow distributions: disentangling the influence of viscosity and initial 
state fluctuations 


With improved detection techniques at RHIC and new data at the higher LHC energies (where each collision 
creates more than twice as many particles as at RHIC), higher order anisotropic flow coefficients of particle momentum 
distributions have now also been precisely measured m- The full anisotropic flow spectrum {vn} provides us with 
enough information to disentangle the extraction of 77/s for the QGP from the dynamical consequences of the 
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FIG. 4. Left Panel: Gharged hadrons anisotropic flow Vn compared with the state-of-the-art IP-Glasma + viscous hydrody¬ 
namic (MUSIG) simulations for Pb+Pb at 2.76 A TeV. Middle and right panels: Event-by-event 772,3 probability distributions 
compared with the same model simulations. The plots are taken from m- 
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geometric shape of the initial-state and its quantum fluctuations that were at the heart of the model uncertainty 
seen in Fig. lnani]. Figure shows recent results from (3+l)-dimensional viscous fluid dynamical simulations 
m performed with the currently most advanced initial state model for heavy-ion collisions, the IP-Glasma model 
m- The left panel illustrates that this approach provides a precise and consistent description of the measured 
charged hadron anisotropic flow coefficients, V 2 to V 4 ^ with (77 /s)qgp = 0.2 [16] (recently refined to {r]/s)QQp = 0.18 
HZ])- The middle and right panels of Fig. show that, even more impressively, the IP-Glasma model initial density 
distributions, after viscous hydrodynamic evolution, yield an almost perfect description of the experimentally unfolded 
Vn probability distributions m- Experimental access to the full spectrum of Vn coefficients and to the probability 
distributions of their event-by-event fluctuations has changed the scope for a precise determination of the QGP 
specific shear viscosity: a phenomenological extraction of {r]/s)QGP with a relative precision of order 5-10% now 
appears within reach. The same data, as well as several additional recently identified flow fluctuation and flow 
correlation observables that focus more specifically on the fluctuating behaviour of the flow angles [ 1 HH 2 I] also 
promise a way to simultaneously constrain the spectrum of the initial-state quantum fluctuations. The procedure 
involves extensive modeling combined with advanced statistical analysis tools [22], mirroring modern studies of the 
cosmic microwave background [23] . 
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FIG. 5. Time line of precise extraction of the specific shear viscosity of QGP. The plot is taken from the “Hot & Dense QCD 
White Paper”, solicited by the NSAC subcommittee on Nuclear Physics funding in the US. Available at http://www.bnl. 
gov/npp/docs/Bass_RHI_WP_final.pdf 


In all, through efforts in the development of sophisticated phenomenological models and improvements in the 
precision of the anisotropic flows measurements, the extraction of the QGP r]/s has become increasingly accurate. 
Fig. [^summarizes a time line of the improvement in constraints placed upon r]/s of the QGP over the past decade. 
We can see that the uncertainty of the value of the QGP r]/s has shrunk dramatically. The convergence is due to 
progress on both the theoretical and experimental sides. At the current stage, theoretical models and experimental 
measurements are both beginning to reach the sensitivity necessary to constrain even the temperature dependence 
of {r]/s){T) laillllS], as well as that of other transport coefficients, such as the bulk viscosity [26] and various 
second-order transport coefficients. 

Outlook and Challenges 

While phenomenological modeling with hydrodynamics has led to an excellent description and to successful 
predictions of most soft (pT< 2 GeV/c) hadronic observables in relativistic heavy-ion collisions, there are some 
“outliers” that the current state-of-the-art theoretical approaches have difficulties explaining. 

In the left panel of Fig.[^ the CMS collaboration measured charged hadron anisotropic flow in 0 - 0 . 2 % ultra¬ 
central Pb+Pb collisions at 2.76 A TeV [27]. The relatively similar size of the measured V 2 and vs coefficients is a 
striking feature which challenges current phenomenological model descriptions [28] (while 62 and 53 are predicted 
to be similar in these collisions, shear viscous effects should suppress vs more strongly than V 2 ). Additionally, in 
the past 3 years, a surprisingly large direct photon elliptic flow has been reported by the PHENIX collaboration 
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FIG. 6. Left Panel: Charged hadrons anisotropic flow Vn compared with viscous hydrodynamic simulations for 0-0.2% ultra¬ 
central Pb-hPb collisions at 2.76 A TeV [27]. Middle panel: Direct photon elliptic flow compared with event-by-event viscous 
hydrodynamic models [38]. Right panel: Charged hadrons elliptic flow in peripheral Pb+Pb collisions and central p+Pb 
collisions at LHC compared to the model calculations [39] . 


at RHIC [29] |3T] and later confirmed by the ALICE collaboration at the LHC [30l [32]. All current theoretical 
models underestimate the measured photon momentum distributions and elliptic flow by factors of 2 to 4 [33H37] (a 
typical model/data comparison is shown in the middle panel of Fig. [^. This major challenge has become known in 
the field as the “direct photon flow puzzle”. Finally, the right panel of Fig. shows that the current state-of-the- 
art phenomenological model with IP-Glasma initial conditions m fails to provide a consistent description of the 
measured elliptic (and triangular, not shown) flow coefficients in p+Pb and Pb+Pb collisions at the LHC [17]. While 
this does not necessarily invalidate the hydrodynamic paradigm, it at least shows that our theoretical understanding 
of the internal structure of the proton and of the fireballs it creates when colliding with a large nucleus is still 
incomplete. 

Resolving these puzzles requires further improvements of the theoretical model from the description of initial state 
fluctuations to more realistic bulk dynamic evolution and to the detailed reconstruction of the particle correlations 
in the final states. 

With higher statistics accumulated in the experiments, the anisotropic flow of rare probes, like high-p^ charged 
hadrons {pr > 10 GeV), direct photons, D mesons, J/^p and etc., are now being or will soon be measured to 
satisfactory precision. Compared to the bulk soft hadrons, the anisotropies of these rare particles can probe earlier 
dynamics of relativistic heavy-ion collisions. In principle, they are more sensitive to the temperature dependence 
of rj/s in the high temperature region. A combined analysis of these rare observables with the anisotropic flow of 
charged hadrons at a variety of collisions energies will help us to further constrain our theoretical models and improve 
the precision of our extraction of the QGP transport coefficients. 
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